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ABSTRACT 

O ' We have developed a model which describes the co-evolution of the mass function of 

dense gravitationally bound cores and of the stellar mass function in a protocluster 
^ . clump. In the model, dense cores are injected, at a uniform rate, at different locations 

' in the clump and evolve under the effect of gas accretion. Gas accretion onto the 

cores follows a time-dependent accretion rate that describes accretion in a turbulent 
CO ' medium. Once the accretion timescales of cores of a given age, of a given mass, and lo- 

I cated at a given distance from the protocluster clumps center exceed their contraction 

, timescales, they are turned into stars. The stellar initial mass function (IMF) is thus 

' built up from successive generations of cores that undergo this accretion-collapse pro- 

cess. We also include the effect of feedback by the newly formed massive stars through 
their stellar winds. A fraction of the wind's energy is assumed to counter gravity and 
OO ' disperse the gas from the protocluster and as a consequence, quench further star for- 

I mation. The latter effect sets the final IMF of the cluster. We apply our model to a 

clump that is expected to resemble the progenitor clump of the Orion Nebula Cluster 
(ONC). The ONC is the only known cluster for which a well determined IMF exists 
for masses ranging from the sub-stellar regime to very massive stars. Our model is 
able to reproduce both the shape and normalization of the ONC's IMF and the mass 
, ^ function of dense submillimeter cores in Orion. The complex features of the ONC's 

^ . present day IMF, namely, a shallow slope in the mass range ~ [0.3 — 2.5] Mq, a steeper 

slope in the mass range ~ [2.5 — 12] Mq, and a nearly flat tail at the high mass end are 
reproduced. The model predicts a 'rapid' star formation process with an age spread for 
the stars of 2.3 x 10^ yr which is consistent with the fact that 80 percent of the ONC's 
stars have ages of < 0.3 Myr. The model also predicts a primordial mass segregation 
with the most massive stars being born in the region between 2 and 4 times the core 
radius of the cluster. In parallel, the model also reproduces, at the time the IMF is set 
and star formation quenched, the mass distribution of dense cores in the Orion star 
forming complex. We study the effects of varying some of the model parameters on 
the resulting IMF and we show that the IMF of stellar clusters is expected to show 
significant variations, provided variations in the clumps and cores physical properties 
exist. 

Key words: galaxies: star clusters - Turbulence - ISM: clouds - open clusters and 
associations 
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1 INTRODUCTION 

1.1 ORIGIN OF THE IMF 

The origin of the stellar initial mass function (IMF) and its 
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ing issues in modern astrophysics. Salpeter (1955) made the 
first observational measurement ol the IMF for field stars 
whose masses are > 0.3 M©. He obtained an IMF that 
follows a power law of the form dN/dM = M* ^-^^ (or 
dN/dlogM = M^^'^'"^^ ). Subsequent work by several au- 
thors showed that the IMF possesses a slope of 4>i ~ —0.1 
at low masses (i.e., M < 0.1 M©) mediated by a plateau- 
hke regime with a shallow slope tl^i ~ —0.3 in the mass 
range 0.1 — 0.5 M© and a steeper slope, xph, at larger masses. 
Whereas Miller & Scalo (1979), Scale (1998), and Chabricr 
(2003) argued for a lognormal form for the IMF, Kroupa 
(2002) proposed that the field stars IMF is well described 
by a three component power law. Over the years, it has been 
argued that the IMF of stellar clusters, at least in the inter- 
mediate to high mass end, is 'universal', that is, i/j^ is given , 
within statistical uncertainties, by the Salpeter value and is 
independent of the environment or of the protocluster cloud 
properties. Elmegreen (2008) calculated the statistical un- 
certainties of the IMF slope for randomly sampled IMFs for 
clusters of different masses. He found that the uncertainties 
on iph are of the order of ±0.15 — 0.2 for clusters masses in 
the range 10^ — 3 x lO'' Mq, respectively. However, devia- 
tions from 'universality' at both the low and high mass ends 
have been reported in many observations. At high masses, 
the IMF is observed to be generally top-heavy in young star- 
burst clusters such as Arches (e.g., Figer et al. 1999; Stolte et 
al. 2005; Kim et al. 2006), NGC 3603 (e.g., Harayama et al. 
2008), and R136 in the central region of 30 Dor (Andersen 
et al. 2009). Espinoza et al. (2009) re-observed the Arches 
cluster and obtained an overall slope of iph = —2.1 ± 0.2 
for the cluster and tph = —1.88 ± 0.2 in its inner annulus 
(within 0.2 pc) whereas Stolte et al. (2005) found a slope 
of tph = —1.26 ± 0.07 in the central region and Stolte et al. 
(2005) and Kim et al. (2006) obtained slopes that are in the 
range of — [1 .7— 1 .9] in the second annulus of cluster (depend- 
ing on whether the entire covered mass range is considered 
or split between intermediate and high mass stars). Even for 
less massive clusters, the universality of the IMF is not well 
established. In fact, determinations of the IMF of many star 
clusters show cluster-to cluster's variations in iph and devi- 
ations from the Salpeter value that are in some cases larger 
than those due to statistical uncertainties (e.g., Massey et 
al. 1995; Okumura et al. 2000; Massey 2003; Leistra et al. 
2005,2006; Sharma et al . 2007). Da Rio et al. (2009) found 
that the IMF of the stellar association LH 95 in the Large 
Magellanic Cloud (LMC) has a slope of iph ~ —2 for stellar 
masses M > 1.1 Mg. Many observational groups have at- 
tempted to use their determinations of the IMFs of various 
Galactic and extragalactic clusters in order to favor one or 
the other of the standard IMFs (i.e., the Chabrier IMF and 
the Kroupa IMF). For example, Liu et al. (2009) recently 
derived the IMF of the NGC 1818 cluster in the LMC and 
argued that the IMF of this cluster is consistant with both 
the Chabrier and Kroupa IMFs. A close inspection of the 
IMF of NGC 1818 (Figures 7 and 8 in their paper) clearly 
shows the best fit Chabrier IMF to the data would predict 
a peak at 0.3 Mq, whereas the data indicates that the peak 
of the distribution is around 0.8 Mq . In fact the data shows 
that there is instead a depression at the predicated peak by 
the Chabrier IMF at 0.3 Mq. As for the Kroupa IMF, it 
systematically overestimates the number of stars of masses 
> 1.5 M© in NGC 1818, underestimates the number of stars 



at the peak, and does not go through the error bars of most 
of the points. 

Another related issue is that massive stars appear to 
be preferentially located in the central parts of the clusters 
(e.g., Hillenbrand & Hartmann 1998; Figer et al. 1999; Siri- 
anni et al . 2002; Stolte et al. 2002; Gouliermis et al. 2004; 
Chen et al. 2007, Sharma et al. 2008). The mass segrega- 
tion of stars in clusters could be due to dynamical two-body 
relaxation cS'ects (e.g., Vesperini & Hcggie 1997; Kroupa 
2001; Mouri & Taniguchi 2002; McMillan et al. 2007) but 
also could have a primordial origin. Dib (2007) and Dib et 
al. (2007a,2008a) showed that the shallower-than Salpeter 
IMF in the central regions of the Arches cluster, and of like- 
wise starburst clusters, can be reproduced by models based 
on the rapid coalescence of dense prestellar cores (PSCs) be- 
fore they collapse and turn into stars. Recently. Chatterjee et 
al. (2009) confirmed the results of Dib et al. (2007a) as they 
find that it is unlikely for the mass segregation in the Arches 
cluster, at its observed current level, and for the cluster age, 
to be the result of dynamical mass segregation starting from 
a standard non-segregated Kroupa IMF. They conclude that 
a shallower IMF in the inner regions of the cluster and which 
would be the imprint of the star formation process is needed 
in order to explain the degree of mass segregation presently 
observed in Arches. 

Undeniably, the origin of the IMF is intimately related 
to the origin and evolution of the gravitationally bound 
molecular cloud cores in which stars form. Several observa- 
tional studies using a variety of techniques and wavelengths 
have shown that the slope of the dense prestellar core mass 
function (PSCMF) is well bracketed by the estimated slopes 
of the IMF (Motte et al. 1998, 2000; Johnstone 2000,2001; 
Kirk et al. 2006, Stanke et al. 2006; Johnstone & Bally 2006; 
Alves et al. 2007). The role played by several physical pro- 
cesses on the origin and evolution of the PSCMF has been 
extensively studied both theoretically and numerically. The 
major physical processes considered are: a) gravitational 
fragmentation (e.g., Zinnecker 1984; Larson 1985; Klessen 
et al. 1998; Klessen & Burkert 2001), turbulent fragmen- 
tation (Elmegreen 1993; Padoan 1995; Padoan et al. 1997; 
Padoan & Nordlund 2002; Padoan et al. 2007), gas accretion 
(Zinnecker 1982; Larson 1992; Bonnell et al. 1997; Klessen 
& Burkert 2000; Bonnell et al. 2001a,b; Basu & Jones 2004; 
Bate & Bonnell 2005; Dobbs et al. 2005; Bonnell & Bate 
2006; Clark & Bonnell 2006; Banerjee et al. 2006; Bonnell 
et al. 2007, Heitsch et al. 2008; Clark et al. 2008; Dib et al. 
2008b, Offner et al. 2008, Myers 2009), cores or star coales- 
cence (Field & Saslaw 1965; Nakano 1966; Silk & Takahashi 
1979; Podsiadlowski & Price 1992; Price & Podsiadlowski 
1995; Bonnell et al. 1998; Bonnell & Bate 2002; Elmegreen & 
Shadmehri 2003; Elmegreen 2004; Shadmehri 2004; Davies 
et al. 2006; McMillan et al. 2007; Kitsionas & Whitworth 
2007; Dib 2007; Dib et al. 2007a,2008a), ejection of predom- 
inantly low mass stars (e.g.. Bate et al. 2002; Goodwin et al. 
2004; Bate & Bonnell 2005), and magnetic fields (Shu 2004, 
Dib et al. 2007b; Nakamura & Li 2008; Price & Bate 2008; 
Kunz & Mouschovias 2009). Accretion onto PSCs which are 
on their way to form stars might also be regulated by the 
cS'ects of feedback from protostellar oufiows and jets (Li & 
Nakamura 2006; Nakamura & Li 2007) and ionization fronts 
and winds from massive stars (Adams & Fatuzzo 1996; Dale 
et al. 2005; Dale & Bonnell 2008; Kevlahan & Pudritz 2009). 
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1.2 EVIDENCE FOR ACCRETION IN HIGH 
MASS STAR FORMING REGIONS ? 

Among the above mentioned processes, core coalescence and 
gas accretion might play a central role in modifying the ini- 
tial PSCMF inherited from the gravo-turbulent fragmenta- 
tion of a protocluster clump. Dib et al. (2007a) have shown 
that an efficient coalescence process of dense PSCs in a pro- 
tocluster clump and their subsequent collapse into stars can 
result in a significantly shallower-than Salpeter IMF in the 
intermediate-to-high mass regime. The process based on the 
coalescence and collapse of dense cores can explain the ob- 
served shallow mass functions of young and massive star- 
burst clusters such as Arches, NGC3603 and RC136 (Kim 
et al. 2006; Harayama et al. 2007; Andersen et al. 2009; Es- 
pinoza et al. 2009). In the case of less centrally condensed 
clusters, the coalescence of their precursor PSCs might be 
less important. However, dense PSCs will continue to ac- 
crete gas from their surrounding medium and continue their 
mass growth. Fig. [T] (top) displays the combined PSCMF 
of dense cores in Orion A North, Orion A South, Orion 
B North and Orion B South based on the compilation of 
cores in Orion by Nutter & Ward Thompson (2007). Note 
that the PSCMF displays a tail-like feature at the high-mass 
end (Another complex PSCMF is the one obtained for for 
the W3 region as observed by Moore et al. 2007). However, 
if the high mass cores in the crowded regions OMCl and 
NGC 2024 are excluded, the distribution of dense cores in 
Orion is well fitted with a universal lognormal distribution 
(Andre et al. 2008). Fig. [1] (bottom) displays the IMF of the 
Orion Nebula Cluster (ONC) as initially derived and plot- 
ted by Hillenbrand (1997) and updated by Hillenbrand in 
2003. Both the PSCMF of Orion and the IMF of the ONC 
show an intriguing flattened tail in the high mass regime. 
Furthermore, The IMF of the ONC displays complex fea- 
tures, among which is a plateau in the mass range 0.45-2.5 
Mq, a nearly universal slope between 2.5 and 10 Mq, and 
another plateau in the high mass regime. Note that a tail at 
the high mass end is not observed in the CMF of low mass 
nearby star forming regions as well as in the IMF of young 
embedded clusters which do not contain massive stars such 
as IC 348, NGC 2362 and Chamaeleon I (e.g., Luhman et al. 
2000; Luhman et al. 2003a,b; Lada & Lada 2003; Muench et 
al. 2003; Luhman 2004,2007). 

Based on the estimate of the lower limit of the ONC 
mass and the assumption that the IMF of the ONC can be 
described by a Kroupa IMF, Pflamm-Altenburg & Kroupa 
(2006) argued that the ONC should harbor about 40 OB 
stars that are heavier than 5 Mq whereas only 10 OB stars 
are observed to be present in the cluster. Using N-Body sim- 
ulations, they showed that it would be possible for a compact 
cluster containing about 40 OB stars to eject 30 of them, 
along with many lower mass stars, within a timescale of ~ 1 
Myr. The OB stars should be detected in the mid-infrared 
by the bow shocks they produce as they evade the cluster. 
Runaway stars that could have been ejected from dense clus- 
ters are observed in star forming regions (e.g., Hoogerwerf et 
al. 2000; Gomez et al. 2005). Gvaramadze & Bomans (2008) 
report the observations of three bow shocks produced by O- 
type stars ejected from the NGC 6611 (M16) cluster. Note 
that Huthoff & Kaper (2002) discussed the fact that the 
detection of bow-shocks due to stellar winds from massive 
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Figure 1. Top: Core mass function of tiie Orion star forming 
region. Tlie data points combine tiie cores from the Orion A Nortii 
and South regions and the Orion B North and South regions (from 
Nutter et al. 2007). The full lines include all cores whereas the 
dashed lines correspond to the population of starless cores (i.e., 
no young stellar objects detected within them in the infra-red). 
Bottom: Initial mass function of the Orion Nebula Cluster as 
plotted by Hillenbrand (1997) and updated by L. Hillenbrand in 
2003. 

stars may depend on the local conditions of the ISM where 
the star is moving. However, as compelling as the scenario 
proposed by Pflamm-Altenburg & Kroupa is, there is no 
observational evidence that about 30 OB stars have been 
ejected from the ONC. Even for low mass stars, there are 
also no observational support for the presence of runaway 
low-mass stars in the ONC (O'Dell et al. 2005) contrary to 
what has been suggested by Poveda et al. (2005). 

In this work, we assess how the mass function of PSCs 
is modified by the effects of time dependent gas accretion. 
As initial conditions for the distributions of PSCs in a pro- 
tocluster clump, we use ones resulting from the local turbu- 
lent fragmentation of the clump. We show that the PSCMF 
evolves quickly under the effect of accretion and develops 
a tail-like feature at the high mass end, comparable to the 
one observed for the mass function of dense cores in star 
forming regions which harbor a population of massive cores 
(i.e., Orion). We also account for the transition from dense 
cores/protostars to stars and discuss the evolution of both 
the CMF and IMF as a function of the model parameters. In 
§[2j we briefly discuss the observed properties of star forming 
protocluster clumps and present our prescription for mod- 
eling them. In § 131 we discuss some of the properties of the 
PSCs and in § [4] we present the local mass distributions of 
PSCs that are formed in the clump uniformly over time. In 
§ [S] we describe our adopted accretion model for accretion 
onto the cores, and in § [6] we present the case of a fidu- 
cial model and its application to the Orion Nebula Cluster. 
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Namely, we study the temporal co-evolution of the PSCMF 
and IMF in the fiducial model and show how the model is 
able to reproduce simultaneously the complex features in 
the PSCMF in Orion and the IMF of the ONC. We also 
address the related issues of age spread of stars and their 
mass segregation. The effects of varying the model main pa- 
rameters on the resulting IMFs are presented in § [7] and in 
§ [9l we conclude. All the variables of the model are listed 
and explained in Tab. [T] including the free parameters that 
are studied in § [T] 
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2 PROTOCLUSTER CLUMPS 
2.1 OBSERVATIONS 

Over the last two decades, several studies using a variety of 
wavelengths and techniques have established that star clus- 
ters form in dense (> 10^ cm~^) clumps embedded in a lower 
density parental molecular cloud (e.g., Lada & Lada 2003; 
Shirley et al. 2003; Minier et al. 2005; Allen et al. 2007 and 
references therein). Saito et al. (2007) recently studied, using 
the C^*0 molecular emission line, a large sample of cluster 
forming clumps whose masses and radii vary between [15- 
1500] Mgand [0.14-0.61] pc, respectively. Fig. [2] displays the 
scaling properties for the star forming clumps observed by 
Saito et al. (2007) (based on a re-interpretation of figures 4 
and 6 in their paper). The mass-size, and velocity dispersion- 
size relations which we adopt to further constrain our models 
and which are obtained by performing least square fits to the 
data points in Fig. [2] (the fits are over-plotted to the data in 
Fig. [2]), are given by 



Mc(M0) 
and 

?;c(km s~ 



-j^g3.62±0.14^2.54±0.2f 



(pc), 



1) ^ ^o0.45±0.08^0.44±0.14(p^)^ 



(1) 



(2) 



2.2 MODELS 

In Dib et al. (2007a), we have adopted a protocluster clump 
model that follows an density profile. Albeit this is a 
fairly good representation for star forming clumps, we al- 
low here for generic variations in the clump density profiles 
by assuming that they can be described by the following 
function: 



Pc(r) = 



PcO 



1 + (r/i?co) 



(3) 



where R^o is the clump's core radius, pcO is the density at 
the center, and h a parameter that accounts for variations in 
the clump radial density profiles in their outer regions. For a 
given mass of the clump, the central density is given by the 
following equation, which can be easily solved numerically^ 

PcO = Tff— — , ,^ ■ . (4) 



il^' 4^rV(l + (r/i?,o)')dr' 



^ An analytical solution 
Mc 

47r(l/3)H3Fi^(3/b,l;l+(3/b);-(flc/-Rco)'' ' 
hypergeometric function. 



to Eq. [4] exists and is given by 
where is a (2,1) order 
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Figure 2. Scaling relations for clumps observed in the C^*0 
molecular emission line, namely the mass-size relation (top), and 
the velocity dispersion-size relation (bottom). These figures are 
based on data presented in figures 4 and 6 in Saito et al. (2007). 
Over plotted are least-square fits to the data. 



where is the radius of the clump. Note that a variety 
of profiles have been derived for star forming clumps rang- 
ing from b ' 1.5 to ' 2.5 (e.g., Motte et al. 1998). The 

temperatures of the cluster forming clumps are observed to 
vary between 15 and 70 K (e.g., Saito et al. 2007). We take a 
conservative value of the temperature of T = 20 K, which is 
probably more representative of their central parts where the 
bulk of the mass is located (e.g., Minier et al. 2005). In order 
to further constrain the models and minimize the number of 
parameters, we relate the size of the protocluster clumps to 
their mass using the mass-size relation presented in § 12.11 In 
our models, the proto-cluster clumps are assumed to be in 
equilibrium. Tan & McKee (2002) and Tan et al. (2006) ar- 
gued that cluster-forming clumps can be in a state of equilib- 
rium for a few to several dynamical crossing times. However, 
as our results below will show, the final IMF of a stellar clus- 
ter in a protocluster clump is formed on a timescale which 
is of the order of 0.2 — 0.3 f//,c , where t//,c is the free-fall 
timescale of the clump. Thus, the equilibrium is required 
to hold for much shorter timescales tan those suggested by 
Tan & McKee (2002). Elmegreen & Shadmehri (2003) and 
Shadmehri (2004) assumed that star forming clumps in a 
molecular cloud are virialized. This might be a plausible hy- 
pothesis if the clumps were indeed the dissipative structures 
of turbulence in the interstellar medium. However, numer- 
ical simulations (e.g., Dib et al. 2007b, Dib & Kim 2007) 
show that clumps and cores in molecular clouds are not in 
virial equilibrium. In the absence of detailed information 
about the velocity dispersion inside the cores in the Saito et 
al. (2007) study, we assume that the clump-clump velocity 
dispersion they derived (i.e., Eq|3]) is also valid on the scale 
of the clumps themselves and of their substructure. 
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3 THE PRESTELLAR CORES MODEL 



Whitworth & Ward-Thompson (2001) have apphed a family 
of Plummer sphere-like models to the contracting prestellar 
dense core L1554. This core is similar to the population of 
gravitationally bound cores that can be found in a clump 
and that are considered in this work. They found a good 
agreement with the observations of L1554 if the density pro- 
file of the core has the following form: 



0.0 



pp{rp) = 



PpO 



[l + {rp/Rpoy 



(5) 



where ppo and Rpo are the central density and core radius 
of the PSC, respectively. Note that the radius of the PSC, 
Rp, depends both on its mass and on its position within the 
MC. The dependence of Rp on r requires that the density 
at the edges of the PSC equals the ambient clump density, 
i.e., pp{Rp) = pc{r). This would result in smaller radii for 
PSCs of a given mass when they are located in their inner 
parts of the cloud. The density contrast between the edge of 
the PSC and its center is given by: 



C(r) ^ = ft£ 

Pc[r) pco 



1 + 



r 
RcO 



(6) 



Depending on its position r in the cloud, the radius 
of the PSC of mass M, Rp, can be calculated as being 
Rp{r,M) = a{r) Rpoir,M), where: 



R 



■pO 



^) = (2^) (-tanK.)l - , (7) 



and with a(r) = (C(r)^^^ — 1)^'^^. With our set of parameters, 
the quantity C^^^ — 1 is always guaranteed to be positive. 
The value Rp{r,M) can be considered as being the radius 
of the PSC at the moment of its formation. The radius of 
the PSC will decrease as time advances due to gravitational 
contraction. Both observational (Lee & Myers 1999; Jessop 
& Ward-Thompson 2000; Kirk et al. 2005; Hatchell et al. 
2007; Ward-Thompson et al. 2007) and numerical (Vazquez- 
Semadeni et al. 2005a; Galvan-Madrid et al. 2007; Dib et 
al. 2008c) estimates of gravitationally bound cores lifetimes 
tend to show that they are of the order of a few times their 
free-fall time, albeit decreasing (but still larger than one free- 
fall time) when cores are defined with increasingly higher 
density tracers/thresholds. Thus, we assume that the PSCs 
contract on a timescale, tcont.p which we take to be a few 
times their free fall timescale tff, and which is parametrized 
by: 



tcont,p{r,M) = u tff{r,M) = u 



32 Gpp{r,M) 



1/2 



(8) 



where G is the gravitational constant, v is a. constant ^ 1 
and pp is the radially averaged density of the PSC of mass 
M, located at position r in the clump, and which is calcu- 
lated as being: 



Pp{r,M) = 



Rp{r,M) 



ppO 



-drr, 



(9) 



Rp{r,M) Jo [l + {rp/Rpo)T 

Thus, the time evolution of the radius of a PSC of mass 
M, located at position r in the cloud is given by a simple 
contraction law: 

Rp{r,M,t) = Rp{r,M,0) e-^'/'-" 



(10) 
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Figure 3. The variation of the peak emission in the 850|^m and 
hence of the peak density (in our paper named Ppo) as a function 
of the mass for a sample of dense cores in the observations of 
Johnstone & Bally (2006). A least square fit to the data points 
yields S850,peak OC PpO OC m0-59±0-05. 



One important issue is the choice of the quantity ppo- 
In Dib et al. (2007a), we have assumed that the dense PSC 
have a constant peak density that is independent of the PSC 
mass and position. In this work, we adopt a slightly more 
realistic approach along the following lines: We first assume 
that the minimum density contrast that should exist be- 
tween the center of the PSC and its edge is of the order 
of the critical Bonnor-Ebert sphere value and that is > 15. 
Secondly, we assume that the density contrast between the 
center and the edge of the PSCs depends of their masses 
following a relation of the type: 



ppO OC iW" 



(11) 



Caselli & Myers (1995) found that massive PSCs in the star 
forming regions L1641 (Orion A) and L1630 (Orion B) are 
denser than lower mass ones. The data points in Fig. [3] dis- 
play the peak dust continuum emission at 850n.m, 5*850, peafe 
of dense cores as a function of their masses as observed by 
Johnstone & Bally (2006) in the Orion B molecular cloud. 
The dust continuum emission is directly proportional to the 
dust density, and thus, to the gas density, assuming the gas- 
to-dust ratio is constant and independent of the density. A 
fit to the observed data points in Fig. [3] yields the following 
relation: 



Ss50,peak OC PpO OC M 



M=0.59±0.05 



(12) 



Thus, the density contrast between the center and the edge 
for a PSC with the minimum mass we are considering, Mmin 
(typically Mmin =0.1 Mq) is 15, whereas for a more massive 
PSC of mass M, the density contrast will be equal to 15 x 
{M /M-min,)^ , with p assumed to be in the range [0 — 0.6] . 



6 S. Dib et al. 



4 THE MASS DISTRIBUTION OF 

GRAVITATIONALLY BOUND PRESTELLAR 
CORES 

As initial conditions for the PSCs mass distribution at differ- 
ent clump radii, we adopt distributions that are the result of 
the gravo-turbulent fragmentation of the clump. As in our 
previous work (Dib et al. 2007a), we use the formulation 
given by Padoan & Nordlund (2002, PN02) in order to cal- 
culate the local distributions of PSCs masses. For standard 
parameters characteristic of star forming regions, the PN02 
formalism gives a PSCMF for PSCs that bears a good resem- 
blance to the observations. Hennebelle & Chabrier (2008) 
proposed a derivation of the PSCMF based on the Press- 
Schechter formalism applied to dense cores. Yet, the differ- 
ences between the two models as far as gravitationally bound 
cores are concerned, are quite marginal. In the following, we 
briefly remind what the ingredients of the PN02 model are. 
The model assumes that the probability distribution func- 
tion of an isothermal, turbulent, compressible gas is well 
described by a lognormal distribution (Vazquez-Semadeni 
1994) and is given by: 



P{ln x)d In X - 



: exp 



1 / In X — In a; 

2 I 



d In X, (13) 



where x is the number density normalized by the average 
number density, x = n/n. The standard deviation of the 
density distribution Od and the mean value In x are functions 
of the local thermal rms Mach number, M and Ina; = —(7^/2 
and (j\ = ln(l -I- M^"f'^). PN02 suggested a value of 7 ~ 0.5, 
whereas Kritsuk et al. (2007) using higher resolution simu- 
lation found that 7 ~ 0.260 ± 0.001. The latter value is the 
one adopted in our models. A second step in this approach 
is to determine the mass distribution of dense cores. PN02 
showed that by making the following assumptions: (a) the 
power spectrum of turbulence is a power law and, (b) the 
typical size of a dense core scales as the thickness of the 
post-shock gas layer, the cores mass spectrum is given by: 



N{M) dlog M oc M'^^'-'^-^^d log M, 



(14) 



where j3 is the exponent of the kinetic energy power spec- 
trum, oc fc"*^, and is related to the exponent a of the size- 
velocity dispersion relation in the cloud with /3 = 2q: -f 1. 
However, Eq. [14] can not be directly used to estimate the 
number of cores that are prone to star formation. It must 
be multiplied by the local distribution of Jeans masses. At 
constant temperature, this distribution can be written as: 



P{Mj) dMj = 



:Af 



j^exp 



In Mj - A 




1 ,0 

Mass (in M^^i^ 



10.0 



100.0 



Figure 4. The relationship between the core mass and the ac- 
cretion rate normalization coefficient. This figure is based on the 
combined data of Figure 3 in Schmeja & Klessen (2004). Mq is 
related to the maximum temporal accretion rate by logMo = 
logMrnax+T ■ Overplotted to the data is a linear fit whose param- 
eters are given in the text. 



where 'erf is the error function. The local normal- 
ization coefficient fair) is obtained by requiring that 
jM„a, N(^r,M) dM = 1 in a sheU of width dr, located at 
distance r from the clump's center. Therefore, the local dis- 
tribution of cores generated in the clump, at an epoch t, 
7V(r, M, t), is obtained by multiplying the local normalized 
function N{r, M) by the local rate of fragmentation such 
that: 



7V(r, M, r) dt ■ 



c{r)pcir) 



dt 



< M > (r) tcont,p{r,M) tff^d 



iV(r,M), (18) 



where dt is the time interval between two consecutive 
epochs, < M > is the average core mass in the local distribu- 
tion and is calculated by < M >= / j^^^""" M N{r, M, 0) dM, 
and Ec is a parameter smaller than unity which describes the 
local mass fraction of gas that is transformed into PSCs per 
free fall time of the protocluster clump, f//,ci (i-e., core for- 
mation efficiency per free fall time). In principle, ec might 
have a radial dependence, but for the sake of simplicity, we 
shall assume to be a constant, independent of radius. In 
our model, since PSCs are generated uniformly over time in 
the clump, tc is also independent of time. 



dMj, (15) 



where Mjo is the Jeans mass at the mean density n. Thus, 
locally, the number of cores is given by: 



N{r,M) dlog M = /o(r) M 
P{Mj)dMj 



-3/(4-/3) 



dlog M, 



(16) 



Eq.[T6]can be solved analytically yielding the following 



form: 



N{r,M)dm = fo{r) 



1 + erf 



AlnM 



2^/2ad 



A/^/(^-«dm,(17) 



5 THE ACCRETION MODEL: ACCRETION IN 
A TURBULENT MEDIUM 

Self-gravitating pretstellar cores can accrete mass from their 
surrounding environment (e.g., Klessen 2001; Padoan et 
al. 2005; Schmeja & Klessen 2004, Dib et al. 2008b). As 
they collapse, they also produce bipolar outflows and jets 
which will tend to reduce the accretion rate by removing 
a fraction of the available mass surrounding them (e.g., 
Tomisaka 2002). Thus, the net effect of gas accretion by 
the core and of gas dispersal by its bipolar outflow is given 
by M = Mace — Mioss, where Mace and Mioss are the mass 
accretion and mass loss rates, respectively. Accretion rates 
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onto PSCs are difficult to measure directly from the obser- 
vations. Observations of the accretion rates are mostly ob- 
tained in the protostellar phase and even then, they are usu- 
ally estimated indirectly from the spectral energy distribu- 
tion or eventually from the properties of associated outflows 
(i.e., a correlation between the accretion rate and outflow 
strength, Hartignan et al. 1995; Bontemps et al. 1996; Wolf- 
Chase et al. 2003). The accretion rates are observed to vary 
as the PSC core evolves. Typical accretion rates for Class 
protostars are found to be in the range 10~^ < Mace/ 
M© yr"^ < 10"* (Hartmann 1998; Narayanan et al. 1998; 
Andre et al. 1999; Ceccarelli et al. 2000; Jayawardhana et 
al. 2001; Di Francesco et al. 2001; Maret et al. 2002; Beuther 
et al. 2002a,b), whereas accretion rates of Class I protostars 
are typically an order of magnitude smaller (Henriksen et al. 
1997; Andre et al. 2000) with values ranging between ~ 10~^ 
M© yr"^ and 5 x 10"*^ M© yr"^ (e.g.. Brown & Chandler 
1999; Greene & Lada 2002; Boogert et al. 2002; Young et 
al. 2003). Massive cores that are likely to form O stars are 
observed to have even higher accretion rates in the range of 
10-3 - 10-2 Mq yr-i (Ho & Young 1996; Zhang & Ho 1997; 
Sandell et al. 2005; Beltran et al. 2006; Garay et al. 2007; 
Zapata et al. 2008). 

In the prestellar core phase, Padoan et al. (2005) mea- 
sured the accretion rate of individual cores in numerical sim- 
ulations of turbulent and self-gravitating molecular clouds. 
They argued that the observed accretion rates of prestel- 
lar cores and protostars can be explained by a volume av- 
eraged Bondi-Hoyle accretion rate Mbh (Bondi & Hoyle 
1944). However, it remains unclear whether the Bondi-Hoyle 
formalism, which describes the accretion onto a point mass 
from a homogeneous gas distribution and with no (or a uni- 
form) velocity field, is suited for the case of gas accretion 
by PSCs in a turbulent molecular clump. In a protocluster 
clump, the gas surrounding the PSCs is highly inhomoge- 
neous, and the velocity field could have a rather complex 
topology (e.g., Dib et al. 2007b; Vazquez-Semadeni et al. 
2008). Krumholz et al. (2006) showed that for an ensem- 
ble of accreting objects from a turbulent medium, some 
of the objects will accrete in a fashion that is closer to 
the Bondi-Hoyle mode if the gas surrounding them has a 
well ordered velocity field and a low vorticity level, whereas 
other objects will accrete according to an accretion rate 
in the vorticity dominated regime with an accretion rate 
Mw (derived by Krumholz et al. 2005, and approximated 
by Krumholz et al. 2006). Krumholz et al. (2006) sug- 
gested that the Bondi-Hoyle accretion rate be replaced by 
Mturb ~ [MB_ff(x)-2 -I- M„(x)-2]-^/2. However, the accre- 
tion rate formula proposed by Krumholz et al. (2006) does 
not take into account the self-gravity of the gas. Klessen 
(2001) derived time dependent accretion rates onto cores 
embedded in a self-gravitating gas and proposed to describe 
them by empirical fit functions with a dichotomy of the fit 
parameters covering four distinct mass ranges. Schmeja & 
Klessen (2004, SK04) measured the time dependent accre- 
tion rates onto cores in their simulations of turbulent and 
self-gravitating molecular clouds. Using data from a large 
ensemble of cores (i.e., ~ 6000 cores), they proposed an em- 
pirical fit to their measurements of the accretion rates which 
has the following functional form: 

\ogMsK04iM,t) = log Mo — t e-'/""", (19) 



where Tacc is a parameter which describes the timescale over 
which the accretion rate declines from its maximum value 
and Mo is a parameter which is linked to the mass of the 
accreting PSC at the end of the accretion process. SK04 
found Tacc to be smaller than the PSCs free fall time tff, 
and by averaging over the populations of cores in different 
mass ranges, they obtained a value of Tacc ~ iff/'i- The sec- 
ond parameter. Mo, is related to the maximum value of the 
accretion rate by log M™^ = log Mo - 7. SK04 did not pro- 
vide a relationship between MJlt^ and the temporal mass of 
the cores. Instead, in figure 3 of their paper, SK04 plotted 
the values of M^u^ as a function of the mass of the core at 
the end of accretion process. However, that for many cores 
what SK04 call Mend is not necessarily the final mass of the 
cores as they can be in a phase of active accretion when 
the simulation was terminated. The data in their figure 3 
corresponds to a sample of molecular cloud simulations in 
which turbulence was driven with a variety of Mach num- 
bers (i.e., between 0.1 and 10) and of driving length scales 
(between half and an eighth of the box size). In Fig. |3]we 
combined all the data points of figure 3 of SK04 and fitted 
the Mmax — Mend relation. We found : 

»> -I n-4.78±0.008 ,^0.65±0.17ni -1 /on^ 

Mmax = 10 Mend M©yr . (20) 

By adopting the time dependent accretion rate formula 
of SK04!I (i.e., Eq[T9}, we make the approximation that for 
an accreting core, the actual mass M{t) at each epoch is a 
new final mass which requires a new normalization of the 
maximum accretion rate and thus of Mq. Finally, note that 
the accretion rates in Schmeja & Klessen (2004) are rela- 
tive to an average density of 10^ cm~'^ which is the average 
number density adopted in their simulations. Thus, in order 
to calculate the accretion rates of PSCs located at different 
positions in a protocluster clump, it is necessary to take into 
account the effect of the varying background density. Since 
accretion is directly proportional to the external density, it 
is necessary to scale the accretion rate of a PSC of a given 
mass located at a position r in the cloud where the local 
number density is n(r) by the value of the accretion rate at 
the 10^ cm-3 number density such that: 



^ Schmeja & Klessen (2004) describe tlie collapsing cores (which 
are replaced by a sink particle) in their simulations as protostel- 
lar wliereas we call our cores at the moment they are formed in 
the clump as prestellar. Note however, that the prestellar cores we 
consider in this work are contracting very rapidly (Eg llOU and will 
go through a protostellar phase before forming stars. The cores 
both in our work and in SK04 are not resolved, and thus, the ac- 
cretions rates they derive for accretion onto the sink particles in 
their simulations can be safely applied to describe accretion onto a 
prestellar cores or onto a protostellar cores (protostar-|-envelope) . 
Whether they are termed prestellar cores, protostellar cores, or 
simply cores is not very relevant. The underlying physics is the 
same in the sense that these are point-like objects that are ac- 
creting gas in a turbulent and self-gravitating medium that can is 
either available in the immediate neighborhoud or that is accreted 
from further away. Would the sink particle in SK04 be resolved 
such as to have a real protostar, a disk, and an envelope, one 
could then consider the accretion separately on the protostar or 
on the protostar -|-envelope, bearing in mind that these two quan- 
tities might simply be directly connected. This is however not the 
case, and the distinction need not to be made. 
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Figure 5. A schematic figure showing the set up used in this 
paper. The circles of different colors represent cores of diff'erent 
ages (i.e., injected at difi^erent epochs) along with their local cu- 
mulative mass function. The mass function of cores in the entire 
proto-cluster clump is the sum of all the local distributions. 



Macc{r = r(n,M,t)) = 



n(cm ^) 
105(cm-») 



MsK04{M,t). 



(21) 



Finally, PSCs also re-inject into the cloud or the larger ISM 
a fraction of their accreted mass in the form of outflows and 
jets (Nakano et al. 1995; Matzner & McKee 2000; Tomisaka 
2002; Machida et al. 2007) . Thus, it is also important 
to include the mass loss by outflows . We moderate the 
accretion rates by introducing a mass loss rate due to 
outflows rate such that AIioss ~ 0.1 Mace, which leads to an 
effective accretion rate Macc.eff ~ 0.9 Mace, where Mace is 
given by Eq. [211 



6 THE CO-EVOLUTION OF THE 

PRESTELLAR CORE MASS FUNCTION 
AND THE IMF 

If the accretion rate Macc.eff is generally a function of time 
and of the mass of the PSC, it can be easily shown that the 
time variation of a population of PSCs of mass Al, located 
at a distance r from the center of the clump, and that was 
injected into the protocluster clump at the epoch r, is given, 
at time t, by: 



f dN{r,M, r, t) 



dt 



f^\Macc^eff^('-^)N 



\dM J 



V dM 



(r,M,T,t). (22) 



Whenever a population of PSCs of a given mass M, 
located at a distance r from the center of the cloud has 
evolved (and accreted) for a time that is equal to its con- 
traction timescale, this population of PSCs is collapsed into 
stars and the accretion process onto these new born stars 
is terminated. In the present model, we do not take into 



account any further potential sub-fragmentation of the col- 
lapsing core. Thus, a single core leads to the formation of 
a single star/star system and as a consequence the derived 
IMFs can be compared to system IMFs rather than to sin- 
gle star IMFs. This approach is similar to the one adopted 
in Dib et al. (2007a) for the study of the effects of cores 
coalescence on the mass spectrum in which the criteria to 
turn a PSC into a star was a comparison of the contrac- 
tion timescale to the time dependent coalescence timescale. 
The difference between the approach adopted in Dib et al. 
(2007a) and the one adopted in the present model is the fol- 
lowing: In Dib et al. (2007a), the local fraction of mass of the 
protocluster clump that was injected into PSCs was a quan- 
tity that was flxed at the beginning of the model and there 
was no new generations of PSCs that were injected beyond 
the initial timestep. In the present model, PSCs with a given 
mass spectrum dictated by the local dynamical conditions 
are uniformly injected over time according to the specified 
rate of PSCs formation per unit free-fall time of the proto- 
cluster clump. Thus, the local number of PSCs of a given 
mass, at a given epoch, is the sum of all the local popula- 
tions of PSCs of the same mass that have been injected at all 
epochs that are anterior or equal to the considered epoch (a 
schematic figure of the clump and the populations of cores 
is shown in Fig.[5|. Note that the local populations of PSCs 
of various ages are evolved separately as they are each in a 
different phase of their accretion history, and that they will 
collapse and form stars at various epochs. Thus, Eq. [22] is 
simultaneously solved for all the populations of PSCs whose 
ages are older or equal to the current epoch. The total local 
number of PSCs of a given mass M, at a time t, will be given 
by: 



N{r, M,t)^J2 ^(^' 



(23) 



As in Dib et al. (2007a), we also assume that only a 
fraction of the mass of a PSC ends up locked in the star. 
This implies that a fraction of the mass of the core is re- 
dispersed into the protocluster clump. In addition to the 
effect of the protostellar outfiows, the rest of the mass is 
lost when the stars are formed by the effects of radiation 
pressure and stellar winds. We account for this mass loss in 
a purely phenomenological way by assuming that the mass 
of a star which is formed out of a PSC of mass M is given by 
M* = ^M, where C s? 1- Matzner & McKee (2000) showed 
that ^ can vary between 0.25 — 0.7 for stars in the mass range 
0.5 — 2 Mq. It is unknown whether this result holds at higher 
masses. However, the similarity between the IMF and the 
dense cores mass function observed by Alves et al. (2007) 
in the Pipe Nebula might be an indication of a constant ^ 
across the mass spectrum (i.e., in their case it is ^ ~ 1/3; 
albeit it should be noted that the cores in the Pipe nebula 
are not believed to be all gravitationally bound). In the ab- 
sence of strong observational and theoretical constraints, we 
shall assume that ^ is independent of the mass. We assume 
that the accretion process onto all PSCs in the protocluster 
clump is terminated whenever the kinetic energy exceeds the 
gravitational energy of the clump, resulting in the gas being 
expelled from the protocluster or at least that star formation 
is severely decelerated. Similarly to Dale & Bonnell (2008), 
we consider that only stars whose mass M* exceeds 10 Mq 
loose mass with a mass loss rate, M*, given by: 
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Figure 6. Time evolution of the pre-stellar core mass function (left), and stellar mass function (right) in the protocluster clump with the 
fiducial model parameters. The stellar mass function is compared to that of the Orion Nebula Cluster (ONC, Hillenbrand 1997; updated 
by Hillenbrand 2003). The combined dense core mass function of the four regions in Orion which contain dense cores is also shown (left 
column). The best fit to the data of the ONC is reproduced at t ~ 0.21 t//,c) where t//,c is the free-fall time of the protocluster clump 
and which corresponds to the epoch at which gas is expelled from the protocluster clump (thick line in lower right panel). The size of the 
mass bin in the model has been re-adjusted such as to resemble the observations mass bins. The PSCMF of the model has been scaled 
up by a factor of 15 to account for the difference between the clump mass considered in the model (10* M©) and the combined mass of 
the regions A and B in Orion which is ~ 1.5 - 2 x 10^ Mq (e.g. Nutter & Ward-Thompson (2007). 
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M* = 10" 



30 Mo 



M 



yr 



(24) 



We also consider that the terminal velocity of the wind 
is given by Vinf = 10^ km s~^. The kinetic energy from winds 
is thus calculated as being 



t'=t 



f N{m)M^,{m)vf^t 



m=lOM0 



dm] dt'. (25) 



We assume that only a fraction of E^i^d will be trans- 
formed into systemic motions that will oppose gravity and 
participate in the evacuation of the bulk of the gas from the 
proto-cluster clump. The rest of the energy is assumed to be 
dissipated or carried away from the protocluster clump by a 
small fraction of the mass, particularly if the massive stars 
are not born in the center of the clump (see confirmation of 
this in ij |6.3|) . Thus the relevant energy, Ek,wind, is given by: 



-Efe ' 



— K Eu 



(26) 



where k is a quantity ^ 1. It is very difficult to estimate k as 
its exact value will vary from system to system depending 
on the number of massive stars, their locations, and the 
interactions of their winds. It is a quantity that can only 
be determined by numerical simulations. As a conservative 
guess for the fiducial model, we take k = 0.1. Ek,wind is 
compared at every timestep to the absolute value of the 
gravitational energy. Eg, which is calculated as being: 



Pc{r)^r'^dr, 



(27) 



where pc is given by Eq. O Note that in this work we only 
take into account feedback from massive stars in the form of 
stellar winds. Another important source of feedback from 
massive stars is their ionizing radiation (e.g., Whitworth 
1979, Dale et al. 2005, Lee & Chen 2007; Gritschneder et 
al. 2009, Bisbas et al. 2009) which we intend to include in a 
future work. Finally, we account for the possible modifica- 
tion of the stellar mass function by the effect of stellar winds 
at the high mass end. The variations in the IMF at the high 
mass end will be given by: 



dN^{r,M,t) 



dt 



(r, M,t), 



(28) 



where M* is the stellar mass loss rate given by Eq. 1191 and 
A'*(r,M,t) is the local number of stars, at time t, of mass Af*. 
The model variables and free parameters are summarized in 
TabU] 



A FIDUCIAL MODEL AND COMPARISON 
TO THE ORION NEBULA CLUSTER 



6.1 



As stated in §[E2l the comparison of our models to the Orion 
star formation region is motivated by the fact that Orion is 
the only nearby star forming region which is harboring mas- 
sive cores, in addition to be relatively well sampled in the 
regime of low mass cores. The other nearby star forming 
regions (e.g., Ophiucus, Perseus, Taurus, Pipe Nebula) and 
for which a PSCMF has been determined, are not known 



to host massive cores that could be the progenitors of mas- 
sive stars. Most importantly, on the stellar side, the IMF of 
the ONC has been obtained by a detailed spectroscopic sur- 
vey (Hillenbrand 1997; Hillenbrand & Carpenter 2000) and 
thus, the features of the IMF of the ONC are much more 
reliable then IMF determinations based on photometric sur- 
veys. The IMF of the ONC also sample a mass range that 
extends from The data of the PSCMF of Orion displayed 
in Fig. [T] (top) is the result of the combination of several 
sub regions of Orion where dense cores are present (Orion 
A North, Orion A South, Orion B North, Orion B South) 
and the local dynamical conditions might be different from 
one region to another. Thus any theoretical model of the co- 
evolution of PSCMF and of the IMF in Orion should focus 
on reproducing, as its primary goal, the IMF of the ONC 
only as it is the best constrained observable quantity. 

In this section, we discuss the time evolution of the 
PSCMF and the transition to the IMF for a fiducial model. 
The parameters for this model are: the mass of the proto- 
cluster clump is Mc — 10* Mq, the temperature of the gas 
is assumed to be T = 20 K, the clump's core radius and ra- 
dial density profile exponent are Rco ~ 0.2 pc and b = —2, 
respectively, the exponent of the peak density-mass relation 
is /X = 0.2, the mass fraction of the clump mass that is trans- 
formed into PSCs per free-fall time is = 0.07, the expo- 
nent of the velocity dispersion-size relation is a = 0.44 (from 
Saito et al. 2007), the ratio of the contraction timescale to 
the free-fall time of the cores is = 1.8. and the fraction 
of the mass of the cores that ends up locked into stars after 
they collapse is ^ = 0.1. A value of ^ = 0.1 is motivated by 
the fact that the peak of the PSCMF function in Orion is 
~ 2 M0, whereas the peak of the IMF of the ONC is lo- 
cated at ~ 0.2 Mq. With a mass of 10'' Mq, the clump has 
a radius Rc = 1.42 pc according to Eq. [T] 

We solve the model's equations using a finite differences 
scheme on a (180,180) linear and logarithmic grid in ra- 
dius and mass, respectively. The timestep is chosen to be 
dt = f//,c/300, where f//,c is the initial free-fall time of the 
protocluster clumjjfl. At any given epoch t — i x dt (where 
i is an integer ^ 300) , we solve i times Eq. 1231 Fig. [S] dis- 
plays the time evolution of the PSCs populations in the en- 
tire protocluster clump (left column) in the fiducial model. 
Once PSCs of a given mass Al, of a given age r, located at 
a given distance r from the center of the clump collapse to 
form stars, they are transferred into the IMF (right column). 
For comparison, we also over-plot to the model the PSCMF 
of the four combined regions of Orion (left column) and the 
IMF of the ONC which we aim to reproduce (right column). 

The overall effect of accretion is to shift the character- 
istic mass to higher masses and also to create more massive 
cores (and as a consequence more massive stars) then were 
present if the core mass function was a simple pile up of cores 
formed at different epochs with no accretion involved. How- 
ever, as cores have finite lifetimes, they will be turned into 
stars as time goes by, thus emptying the corresponding bins 
in the PSCMF. For a given population of cores born at the 
same time, and depending on the value of fi, more massive 



^ We have checked, for our adopted timestep size, that a grid size 
of > 120 cells in mass is needed to ensure the convergence of the 



results 
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Figure 7. Time evolution of the ratio of tlie kinetic energy gen- 
erated by stellar winds and the gravitational energy of the proto- 
cluster clump. Time is shown in units of the protocluster clump 
free fall timescale t//,c- The horizontal dashed corresponds to 

Ek,vjind/Egrav = 1 with K = 0.1. 



cores will collapse faster than their lower mass counterparts 
(for fj, > 0). This will lead to a concentration 

In the early phcises, the effects of accretion are not vis- 
ible on the PSCs mass function. The only effects at this 
stage are a pile-up of the consecutives PSCs populations 
that are being injected uniformly as time advances. When 
t ~ 0.04 — 0.06 t//,c, the effects of accretion become visible 
with the widening of the plateau around the peak value and 
the generation of a larger fraction of more massive cores. By 
t ~ 0.08 — 0.11 tff^c, the plateau has advanced until masses 
of Af ~ 20 - 25 M0. A peak in the PSCs distribution is 
also formed at t ~ 0.12 - 0.15 t//,c at M ~ 0.5 Mq. At 
t ~ 0.08 t//,c, the first generations of PSCs collapse into 
stars and the IMF becomes populated (right column) . Since 
in this fiducial model, the more massive PSCs are more cen- 
trally condensed than their less massive counterparts (i.e., 
/i = 0.2 > 0), this leads to the formation of the intermedi- 
ate mass stars first from the first generations of cores before 
low mass stars are formed (there are at this stage no mas- 
sive cores to form massive stars yet). Accretion continues 
to affect the subsequent generations of intermediate mass 
cores leading to the formation of massive cores and to the 
flattening of the PSCs mass function at the high mass end. 
As time advances, the PSCs mass function continues to flat- 
ten at the high mass end and develops a tail-like structure 
for PSCs masses > 50 Mq. In the lowest left quadrant, the 
PSCMF of Orion has been scaled down by a factor of 15 
which is roughly the scaling factor between the mass of the 
protocluster clump we consider here (10* Mq) and the com- 
bined masses of the four Orion regions (1 — 2 x 10^ Mq, 
e.g., Hillenbrand 1997). Although, we have argued earlier 
that an exact match between the combined PSCMF of the 
four Orion regions and of the models is not to be expected 
because of variations that may exist in the local distribu- 
tions in each of those Orion regions, the agreement between 
the observations and the models is surprisingly very good at 
t ~ 0.21 tff,c. 



Figure 8. The IMF at t = 0.21 iff^^m three different regions 
of the protocluster clump in the fiducial model: In the central 
region between the center and the core radius (full line) , between 
one time and two times the core radius (dotted line), between 
two times and four times the core radius (dash-dotted line), and 
in fore regions outside four times the core regions (dashed-triple 
dot line; there is is no stars however in those outer regions). The 
core radius in this model is Rco = 0.2 pc. Note that there are 
no stars present beyond 4 RcO- The adopted bin size in mass is 2 
times bigger than the adopted mass bin in the calculation. This is 
done to enable a better viewing of the global features and smooth 
any sharp bin-to-bin variations. The same bin size is adopted in 
FiES. [T0l[T2lll4ll5l and [13 



At the same time, as all the bins of the IMF are popu- 
lated, at t ~ 0.21 t//,c, the model IMF reproduces almost all 
of the features of the present day IMF of the ONC. Namely, 
a shallow slope in the mass range ~ [0.3 — 2.5] Mq, a steeper 
slope in the mass range ~ [2.5 — 12] Mq, and nearly flat tail 
at the high mass end. The PSCMF in the model and subse- 
quently the IMF would continue to evolve as time advances. 
However, at ~ 0.21 t//,c, the kinetic energy generated by 
stellar winds becomes comparable or larger than the binding 
gravitational energy in the protocluster clump (see Fig. [T]). 
This would inevitably lead to the ejection of the gas from 
the protocluster clump and to the settling of the IMF into 
a form that is very similar to the present day mass function 
of the ONC. The dispersal of the gas being dispersed from 
the central regions is observed with a layer of star forming 
molecular gas (Genzel & Stutzki 1989), standing between 
an outwardly moving ionized gas produced by the central 
massive stars (O'Dell et al. 1994) and thin layer of neutral 
gas (O'Dell et al. 1992). 

An important aspect of our model is that it reproduces 
simultaneously both the shape and the normalization of the 
present day IMF of the ONC and the mass function of dense 
sub-millimeter cores in Orion. Note that the synchronization 
that occurs between the setting of the final IMF and of the 
gas ejection from the protocluster clumps at t ~ 0.21t//^c is 
related to our choice of k = 0.1. However, if k had a differ- 
ent value in the range 0.1 — 0.5, the epoch of gas dispersal 
from the cluster would shift to t ~ 0.18 — 0.2 t//,c, when 
the model IMF is already in good agreement with that of 
the ONC. The most massive star in the ONC of mass ~ 50 
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Mq is not reproduced by the model. This could simply due 
to very specific conditions in the region of the protocluster 
clump where this star is born and that are not taken into 
account by the model. A process that could produce a very 
massive core from which this massive O star could be born is 
the coalescence of two or more low mass cores as suggested 
by Dib et al. (2007a) (see also Edgar & Clarke 2004). It 
could otherwise be due to the inadequacy of the accretion 
rate normalization that we have used in Eq. [2T]to describe 
the accretion process at the very high mass regime, where 
the effects of radiation from the nascent protostars might 
be important in regulating the accretion rate from the sur- 
rounding medium. Overall, this fiducial model based on the 
co-evolution of the PSCMF and of the IMF is able to re- 
produce most of the observed features of both the PSCMF 
in Orion and the IMF of the ONC. Furthermore, the model 
suggests that the present day mass function of the ONC 
has a primordial origin. Thus, the results of our model ar- 
gue against the dynamical scenario proposed by Pflamm- 
Altenburg & Kroupa (2006). 



6.2 AGE SPREAD OF STARS 

Our proposed scenario for star formation in the ONC based 
on the co-evolution of the PSCMF and of the IMF through 
gas accretion by the cores, their collapse to form stars and 
the quenching effect of star formation by stellar winds im- 
plies a very small age spread for stars in the ONC. In the 
fiducial model the first stars are formed at t ~ 0.07t//,c and 
the star formation is quenched at t ^ 0.21f//,c. The age 
spread in the fiducial model is of the order of ~ 0.14 i//,c ~ 
0.14 X 1.62 X 10*^ yr ~ 2.3 x 10^ yr, where ~ 1.62 x'lO® 
yr is the free-fall time of the protocluster clump with the 
assumed model parameters. In the ONC, the measured age 
spread is ~ 3 Myr based on age determinations made by 
fitting isochrones to the ONC's color magnitude diagram 
for stars with masses < 4 Mq (Hillenbrand 1997). Palla & 
Stabler (1999) argued that the ONC hosts stars that are as 
old as 10 Myr with the bulk of the stars having an age of ~ 2 
Myr. However, Hartmann (2003) pointed out that the de- 
termination of stellar ages using masses in the range 0.4 — 6 
Mq, which was also used by Tan el al. (2006) to argue for 
star formation occurring over several crossing times in the 
ONC, tend to systematically overestimate the ages of the 
stars because the birth line age corrections have been under- 
estimated. However, as pointed out by Hillenbrand (1997), 
about 80 percent of the stars in the ONC have an age that 
is < 0.3 Myr. Furthermore, based on their observations of 
high correlations between the motions of the residual gas 
and of stars in the ONC, Fiiresz et al. (2008) argued that 
the entire system must be very young with an age of at most 
one crossing time. Thus, our model reproduces the charac- 
teristic age spread of most of the stars in the ONC. The 
other 20 percent of stars in the ONC with estimated ages 
> 0.3 Myr may possibly have their ages overestimated, be 
foreground stars (Hartmann 2003), or eventually may have 
formed during the early assembly phase of the protocluster 
clump in a scenario similar to the one proposed by Burkert 
& Hartmann (2004) and Hartmann & Burkert (2007). 
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Figure 9. Time evolution of the ratio of the wind energy to the 
gravitational energy in models that have different values of the 
exponent of the peak density-mass relation fj.. The value of = 0.2 
(black line) corresponds to the fiducial model. The vertical line at 
t ~ 0.2ltff i. marks the epoch at which the IMFs of the different 
models are compared. The horizontal dashed line corresponds to 

Ek,mind/Egrav = 1 with K = 0.1. 
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Figure 10. The IMF at t = 0.21tff ,. in models which have 
different values of the exponent of the peak density-mass relation, 
fi. The value oi /i = 0.2 (black line) corresponds to the fiducial 
model. 

6.3 MASS SEGREGATION 

Observations of stellar clusters often report evidence for 
mass segregation with the most massive stars being pref- 
erentially located in the inner regions of the cluster (e.g., 
Pandey et al. 1992; Subramanian et al. 1993; Malumuth & 
Heap 1994; Brandl et al. 1996; Hillenbrand & Hartmann 
1998; Fisher et al. 1998; Figer et al. 1999; Sagar et al. 2001; 
Stohe et al. 2002; Le Duigou & Knodelseder 2002; Sirianni et 
al. 2002; Lyo et al. 2004; Gouhermis et al. 2004; Sharma et al. 
2007,2008). Using N-body simulation with mass-segregated 
and non-mass segregated clusters, Baumgardt et al. (2008) 
showed that primordial mass segregation explains better the 
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correlation between the slope of the stellar mass function 
and the clusters degree of concentration as observed by De 
Marchi et al. (2007). In the particular case of the ONC, an 
inspection of the spatial distribution of stars in the cluster 
(Fig. 3 in Hillenbrand et al. 2007), shows that two out of 
the six most massive stars are located, in projection, at a 
distance r > 0.5 pc from the cluster center, and one star 
is located, in projection, at a distance of > 1 pc from the 
center, whereas the three other stars are located, in projec- 
tion, within a region r < 0.3 pc from the cluster's center. 
Bases on time scale estimates Bonnell & Davies (1998) ar- 
gued that the massive stars in the ONC must have formed 
in or near the centre and were not the result of dynamical 
mass segregation. 

As mentioned in §[0] the mass segregation of stars in a 
cluster can have a primordial or a dynamical origin, or both. 
In this section, we quantify the effects of primordial mass 
segregation in our accretion-collapse-feedback (ACF) model 
using the above described fiducial model. Fig. [8] displays the 
IMF of the fiducial model in different regions of the cluster, 
namely, in the inner region within one core radius (i?cO = 0.2 
pc), in the second annulus between one and two times the 
core radius, between two and four times the core radius, and 
for the outer region. Fig. [8] shows that albeit the bulk of the 
stars with masses < 20 M© are found in the inner region 
(i.e., r < Rco), the slope of the IMF flattens in the interme- 
diate to high mass regime when going from the inner to the 
outer regions. The most massive stars are born in the region 
located between [2 Rco — 4 i?co] = [0.4 — 0.8] pc. The location 
of the massive stars in the model is in good agreement with 
the positions of at least three of the six most massive stars in 
the ONC. Aside from the fact that one or more of the three 
stars which are within a projected distance of < 0.3 pc from 
the ONC's center might be in reality at a larger physical dis- 
tance from the center in the three-dimensional space, there 
are also other ways to account for the apparent existence of 
these three massive stars in the center of the ONC. One pos- 
sibility is that those three stars have formed, like the other 
ones, in the region 2 Rco — 4 i?cO and sinked to the center 
of the cluster by dynamical interactions (e.g., Allison et al. 
2009). The second possibility, following ideas developed in 
Dib et al. (2007a), is that those few stars were born from the 
coalescence of cores in the center of the proto-ONC clump 
where the latter are expected to be closely packed and where 
the coalescence of cores is expected to be more efficient. The 
coalescence of cores will allow the formation of more massive 
ones and subsequently, after these cores have collapsed, the 
formation of the massive stars directly in the central region. 
The reason why the most massive stars are formed in the 
region [2 Rca — 4 i?co] = [0.4-0.8] pc in our model and not in 
the very inner region lies in the fact that for massive cores 
of equal masses, the cores located in the very inner region 
are smaller, as discussed in § [31 and thus they have larger 
average densities and shorter contraction timescales. The 
shorter contraction timescales of cores in the inner regions 
will lead them to collapse fatser and thus, have a shorter 
accretion history which will hinder their development into 
more massive cores and subsequently more massive stars. 
Overall, the mild mass segregation observed in the ONC 
and in our model is most probably another indication of the 
young age of the cluster. 
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Figure 11. Time evolution of the ratio of the wind energy to the 
gravitational energy in models that have different values of the 
ratio of the contraction timescale to the free fall time scale of the 
cores V = tcont/tff- The value of i/ = 1.8 (black line) corresponds 
to the fiducial model. The vertical line at t ~ 0.21 t//,c marks the 
epoch at which the IMFs of the different models are compared. 
The horizontal dashed corresponds to Ef^ wind/ Egrav = 1 with 

K = 0.1. 
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Figure 12. The IMF at t = 0.21 t//,c in models where the ratio 
of the contraction timescale to the free fall time scale of the cores 
^ = icont,p/'tff = 1.2,1.5, and 1.8 and at t = 0.17 t//,c foi* 
n = 2.4 and t = 0.19 t//,^ for n = 3. The value of i/ = 1.8 (black 
line) corresponds to the fiducial model. 



7 PARAMETER STUDY 

In this section, we investigate the effects on varying some of 
the model's main parameters that are listed in Tab.[l]on the 
resulting IMF. However, for the reader conveniency, we first 
summarize in S 17. li the parameters of the model and explain 
their related expected effects. 
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7.1 Summary of the model parameters 

In this section, we summarize the free parameters in the 
model and their range of variations as constrained by obser- 
vations or theory. We also study the effect of parameter vari- 
ations on the results of the model. The scheme we adopt to 
specify the time dependent accretion rate is not affected by 
free parameters, because it is calibrated using the numerical 
simulations of SK04. The model contains 8 free parameters 
listed in Tab. [1] They can be divided into three categories 
related to the clump properties, the cores properties, or the 
feedback model. Of the 8 free parameters in Tab.[T] the four 
most important ones are labeled IP, while the other four are 
labeled P. 

In the category of clump parameters, the first one is the 
clump mass, Mc. In the fiducial model and in all models be- 
low, Mc has been assigned the same value of Mc — W^ttiq, a 
characteristic mass for protocluster clumps. For a given core 
formation efficiency, changing Mc will in principle only affect 
the vertical normalization of the PSCMF and of the result- 
ing stellar IMF. However, due to the finite size of the clump 
and to the Afc — Rc relation, a more massive clump would 
have a larger radius, which can introduce some minor effects 
due to star formation in the outer regions. The second clump 
parameter is the clump core radius, Rco- In all models we 
chose Rco = 0.2 pc, which is a characteristic value for stellar 
clusters and for their protocluster clump progenitors (e.g., 
Hillenbrand & Hartmann 1998). The third clump parameter 
is the exponent of its density profile in the outer regions, h. 
This parameter can be considered to describe the effect of 
the environment in setting the clump density profile. The 
most quoted value of h in the literature is & ~ 2. However, 
some observations suggest that h may vary from clump to 
clump (e.g. Motte et al. 1998). Thus, in the parameter study 
below, we allow 6 to vary in the range 1.4 — 2.2, which affects 
the accretion rate that is directly proportional to the back- 
ground density. The last of the clump parameters is the rate 
of core formation efficiency per free fall time, t^- Its value 
as a function of Mc is not well constrained by numerical 
simulations or theoretical considerations and may vary from 
clump to clump. Depending on physical conditions, such as 
the magnetic field strength in the clump, the value of tc may 
range from 0.01 or less, to a few times 0.1 (e.g., Nakamura 
& Li 2008; Vazquez-Semadeni et al. 2005b; Padoan & Nord- 
lund 2009). In all the models we choose tc = 0.07. Changing 
Ec would, to first order, only change the normalization of 
the PSCMF and of the IMF and not their shapes. This is 
true as long as the effect of feedback is not considered. If 
feedback is taken into account, then models with larger ec 
values would form a larger number of stars with masses > 10 
Mq (because of the larger normalization in the total mass 
of stars) and thus would have a shorter evacuation time of 
the gas from the protocluster because the wind energy is 
proportional to M*. 

The two free parameters related to the properties of 
PSCs are the contraction timescale of the PSCs normalized 
to their free fall time, v = tcont,p/'tf f , and the exponent 
of the peak density (or concentration)-mass relation of the 
PSCs, /i. As will be illustrated below I7.3p . larger values 
of V imply longer accretion timescales and hence a smaller 
number of low mass cores. This also implies a relatively 
larger fraction of intermediate and high mass cores and thus 



larger fractions of intermediate and high mass stars. As will 
be shown below in § 17.21 the effect of increasing the value 
of fj, is to shorten the free-fall time and therefore the accre- 
tion timescale. Thus for larger ^ values, PSCs of larger mass 
will have less time to accrete compared to lower mass PSCs. 
This results in a decreased fraction of massive stars and a 
steeper IMF at the high-mass end. 

The last two parameters of the model are those related 
to the feedback scheme. The mass fraction of the PSC that 
ends up into a star (the star formation efficiency of an in- 
dividual PSC), ^, does not affect the shape of the resulting 
IMF, nor its vertical normalization; it only shifts the IMF 
horizontally (towards larger or smaller masses). We have as- 
sumed ^ = 0.1 in all models. The parameter k gives the frac- 
tion of the winds energy that opposes gravity in the proto- 
cluster clump. The value of n would not affect the evolution 
of the PSCMF and of the IMF if the stellar feedback were 
not taken into account. The value of n is uncertain and may 
depend on the other parameters. All we know is that k < 1. 
The best fit to the ONC data was obtained for k, — 0.1. 

The radius of the clump, Rc, and the level of turbu- 
lence in the clump described through the exponent of the 
velocity dispersion-size relation in the clump, a, are not free 
parameters because they are constrained by the observed 
Mc — Rc and Rc — Vc relations (Saito et al. 2007; Eq. [T] 
and Eq. [2]). Nevertheless, in 17.51 we varied the value of q, 
within a range allowed by the observational uncertainties in 
those relations, in order to explore the effects of the level of 
turbulence on the IMF. Higher values of a cause a shift in 
the position of the characteristic mass in the PN02 model 
that we adopt as the initial distribution of core masses. 

As the IMF of a stellar cluster is a time evolving dis- 
tribution until star formation is quenched due to the evac- 
uation of the gas from the protocluster clump, ideally, in 
investigating the effects of varying the model's main param- 
eters, we would need to compare the IMF of the different 
models at the time t = texpuise, which corresponds to the 
epoch at which Ek,wind/ Egrav ~ 1, when the gas would 
have been expelled from the protocluster clump. However, 
as stated above, k is uncertain and may vary from system 
to system. If a protocluster clump formed massive stars at 
a slower pace, the condition Ek,wind/ Egrav ~ 1 would be 
reached at a later epoch, leaving more time for other feed- 
back processes not considered in this work, such as jets and 
outflows, to operate. Furthermore, if massive stars appeared 
at later stages, our assumption of clump equilibrium may 
not be justified, and one should perhaps consider the effect 
of clump contraction. To avoid such complications, we shall 
compare the IMFs of the different models at the same epoch, 
independent of whether the gas has been expelled from the 
protocluster clump, totally or partially. With respect to the 
values of the parameters in the fiducial model, we vary the 
values of 6, p, v, and a. 

7.2 Effect of the exponent of the 

mass-concentration relation of the cores 

Fig. [9] displays the time evolution of the ratio of 
Ek^wind/ Egrav iu fivc models where fj, has the values of 0, 
0.1, 0.2, 0.4, and 0.6. The model corresponding to the case 
with the fiducial value {fi — 0.2) is shown with the black 
line. The figure shows that the expulsion of the gas occurs 



Origin of the IMF of stellar clusters 15 



at different epochs in the evolution of the protocluster clump 
for a fixed value of k. Fig. [TO] displays the IMFs for the dif- 
ferent models at t — 0.21 t//,c. The IMF are very similar at 
the low mass end (< 0.1 M©) except for a dip in the IMF 
at ~ 0.05 M0 for the low ^ values. The peak of the IMF is 
located at ~ 0.1 M© for the low /i models and at ~ 0.8 — 1 
M0 in the models with the higher /i values. The models are 
also different in the intermediate to high mass end. Low ^ 
values favor the formation of massive stars and a shallow 
slope in that mass range, whereas the value of is much 
higher in models with higher ^ values. In fact, large values 
of /i imply that the cores are more centrally condensed and 
have smaller contraction timescales than their counterparts 
in models with lower /i values. Thus the intermediate and 
massive cores in the higher /i models have shorter accretion 
timescales than their counterparts with lower /i values and 
therefore proceed faster to form stars before more massive 
cores can be built up by accretion. Thus, larger ^ values, al- 
though, they allow massive cores to collapse faster and form 
a small fraction of massive stars, they delay the formation 
of a sizable fraction of massive stars and consequently delay 
the expulsion of the gas from the protocluster clump as can 
also be seen in Fig. [9] 

7.3 Effect of the contraction timescale of the cores 

Fig. [TT] displays the time evolution of the ratio of kinetic 
energy from stellar winds to the gravitational energy (for 
K = 0.1) for 5 simulations where the ratio of the contraction 
timescale of the cores to their free-fall time is, at all time 
during their contraction, maintained at v — 1.2, 1.8, 2.4, and 
3. Because in the models where the cores contract slower 
(i.e., higher f values, v — 2.4, 3), the feedback becomes very 
important (i.e., Eh,wind/ Egrav ~ 100 for k = 0.1), although 
massive stars appear slightly later than in models with lower 
v values, the gas is expected to be expelled from the proto- 
cluster at an earlier epoch in those models. Thus the IMFs 
for the different models are compared at t = 0.21 t//,c for 
the cases where u — 1.2, 1.4, and 1.8 and ai t — 0.17 t//,c 
for the case with v — 2.4 and at t — 0.19 t//,c for the model 
with u = 3. The IMFs of the models with various u values 
are displayed in Fig. [T^l The IMFs of models with values of 
f higher than the fiducial value show the presence of mas- 
sive stars of increasing mass with increasing f and a clear 
absence of low mass stars. On the other hand, the IMFs 
of models with values of v that are lower than the fiducial 
value show a marked absence of massive stars and similar 
populations of low mass and a decreasing intermediate mass 
population with decreasing v. The accretion history of PSCs 
with lower f values is shorter and inhibits the development 
of a substantial population of massive cores that would pro- 
ceed to form massive stars, whereas cores with high v values 
accrete over more extended timescales and evolves towards 
being more massive with increasing v values, thus keeping 
only a reduced population of low mass core at any given 
epoch and consequently a small fraction of low mass stars. 

Except for the presence of massive stars, the particular 
shape of the IMFs in models with high v values reminds 
of the peculiar IMF in Taurus (Luhman et al. 2003b) (the 
massive stars can disappear from the model IMF if a dif- 
ferent normalization is used, i.e., a lower mass clump). The 
latter exhibits a peak at ~ 0.8 — 0.9 Mq, and a rapid decline 




0.0 0.1 0,2 0.3 0.4 

Figure 13. Time evolution of the ratio of the wind energy to the 
gravitational energy in models that have different values of the 
exponent of the power law profile of the protocluster clumps, b. 
The value of 6 = 2 (black line) corresponds to the fiducial model. 
The vertical line corresponds to the epoch at which t ~ 0.21 iff ^.. 
Note that since the clumps in these models have a different density 
profile, they consequently have different free-fall times t//,c. The 
horizontal dashed corresponds to Sfc^mind/^grai) = 1 with re = 
0.1. 

when going to lower mass stars. Heyer et al. (2008) showed 
that the magnetic field in the Taurus molecular cloud com- 
plex is relatively strong enough and may cause the cloud 
to be, at least for the diffuse part, magnetically subcritical. 
Dense cores that would form in Taurus will thus be likely 
characterized by longer contraction timescales than in other 
regions and thus will be able to accrete over more significant 
periods of time in a similar fashion to the cores with high u 
values in our models. 

7.4 Effect of the protocluster clump density 
profile 

Fig- [131 displays the time evolution of Ek,mind/ Egrav in a se- 
ries of models where a clump of the same mass as the fiducial 
model (i.e., Mc = 10* Mq) possesses a different density pro- 
file characterized by the value of b in Eq. O The values of 
b that are considered are b = 1.4, 1.6, 1.8,2 (fiducial), and 
2.2. Thus, for the same value of Mc, and up to the same 
radius Rc, the clumps with higher values of b will be more 
centrally peaked. The IMFs of the models with different val- 
ues of b are displayed in Fig. [14] at t = 0.21 t//,c and in 
Fig. [TJ] at t = tcxpuise which corresponds to the epoch at 
which Eh,wind/ Egrav ~ 1 (with K=0.1 in all models). 

Whereas the IMF for models for with b — 2, and 2.2 
display a relatively similar IMF, the IMFs of models with 
fe = 1.4 and 6 = 1.6 deviate substantially from the IMF of 
the fiducial model, while the model with b — 1.8 is inter- 
mediate between the two classes. Aside from the issue of 
normalization (the same core formation efficiency of ec = 7 
percent used in all models leads to larger number of cores 
formed in the more centrally peaked clumps), the slope of 
the IMF at the high mass end in the models with 6 = 1.4 
and b = 1.6 are shallower than in the models with higher 
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Figure 14. The IMF at t = 0.21 t//,c in models which have 
different values of the exponent of the power law profile of the 
protocluster clump, b. The value of fe = 2 (black line) corresponds 
to the fiducial model. 



b values (alternatively one can say that the plateau in the 
intermediate mass regime in the low b models is wider). This 
is due to the fact that in all these models, like in the fidu- 
cial one, intermediate mass cores proceed to collapse at a 
faster pace than their low mass counterparts (because A* > 0, 
see Fig. [6]). The intermediate mass cores are thus prevented 
from accreting additional mass and widening the plateau in 
the mass range 0.1 — 2 M©; whereas lower mass core con- 
tinue to accrete and become proportionally more numerous 
in the large b models than in the low b models. When these 
cores finally proceed to collapse, the signature of a steeper 
PSCMF remains imprinted in the generated IMF. If one now 
compares the liMFs of the models with different b values at 
t = texpuise, the fact that fewer massive stars have formed 
in the models with low b values delays the evacuation of the 
gas from the protocluster in those models. Thus, in models 
with lower 6 values, additional generations of low and inter- 
mediate mass cores have time to form and accrete and cause 
the widening of the intermediate mass plateau as compared 
to the case of the high b models. 



7.5 Effect of the turbulence level in the 
protocluster clump 

In the previous sections, we have used the velocity 
dispersion-size relation inside the clumps derived by Saito 
et al (2007). Here we evaluate the effects on the IMFs of 
adopting a different velocity dispersion-size relation expo- 
nent which takes into account the scatter around the cen- 
tral value in the Saito et al. data. Since the value derived 
by Saito et al. (2007) is q = 0.44 ± 0.14, we compare the 
IMF of the model with the fiducial value, a = 0.44 to two 
other models where a — 0.3 and a = 0.58. Fig. 1161 displays 



the time evolution of the ratio Ei. 



d/Eg 



in the models. 



Overall, the evolution of this ratio is very similar in the three 
models, albeit feedback seems to be slightly more important 
in models with a smaller a value. The corresponding IMFs 
are displayed in Fig. [17] at t = 0.21 t//,c for models with 



Figure 15. The IMF at t = iexpeii which corresponds to the 
epoch at which i?fc_u,i„£j/Bg^^^, r^i (for n = 0.1) in models which 
have different values of the exponent of the power law profile of the 
protocluster clump, b. The value of fe = 2 (black line) corresponds 
to the fiducial model. 



& 1 




0.5 



Figure 16. Time evolution of the ratio of the wind energy to the 
gravitational energy in models that have different exponents of 
the velocity dispersion size relation in the clump, a. The value 
of Q = 0.44 (black line) corresponds to the fiducial model. The 
vertical line corresponds to the epoch at which t ~ 0.21ty^ c- 
The horizontal dashed corresponds to E]^ wind/ ^grav = 1 with 

K = 0.1. 



a = 0.44, and a = 0.58, and at t = 0.19 i//,c for the model 
with a — 0.30 since for this model the ratio Ek,wind/ Egrav 
is closer to unity at this epoch. The shapes of the IMFs with 
different values of a are surprisingly similar at the high mass 
end. In the absence of time dependent processes such as ac- 
cretion (which is considered here) , the expected slopes of the 
IMFs would be similar to the slopes of the PSCMF and are 
expected to have the values of ~ —2.25, —2.41, and —2.63 
for a = 0.3,0.44, and 0.58, respectively (PN02). The only 
noticeable difference between these three IMFs and which 
is due to the effect of varying the velocity dispersion-size 
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Figure 17. The IMF at t = 0.21 tff^^ for the models with the 
exponent of the velocity dispersion-size relation a = 0.44 (fiducial 
model) and a = 0.58 and at t = 0.19 f//.c for o: = 0.30. 



relation exponent lies at the low mass end. Higher values of 
Q imply higher Mach numbers and a shift of the character- 
istic mass in the injected PSCs populations towards smaller 
masses (PN02) which is the eflfect that can be observed in 
Fig. [13 



8 COMPARISON TO PREVIOUS WORK 

Several authors have investigated the effects of a number of 
physical processes on the mass function of dense cores or 
on the IMF (see Bonnell et al. 2007 and references therein). 
Among those, analytical derivations of the effects of gas ac- 
cretion on the dense core mass function have been presented 
by Bonnell et al. (2001), Basu & Jones (2004), Bate & Bon- 
nell (2005), and more recently by Myers (2009). Bonnell et 
al. (2001b) presented a two stage model for the evolution 
of the dense core mass function in a protocluster clump. In 
the Bonnell et al. (2001b) model, in the regime where the 
gravitational potential of the cluster is dominated by the 
clump gaseous component, gas is accreted onto a core/star 
whenever it is located at a distance from it that is smaller 
to the tidal lobe-radius. As shown by Bonnell et al. (2001b), 
the result of gas accretion in this regime would be to gener- 
ate a mass spectrum with a slope of tph = —1.5. In a second 
regime where the gravitational potential of the protocluster 
is dominated by stars, the cores/stars would accrete follow- 
ing a Bondi-Hoyle like accretion scheme with the accretion 
radius being equal to the Bondi-Hoyle accretion radius. The 
result of this accretion is to steepen the slope of the result- 
ing IMF inherited from the previous phase to values in the 
range of ~ —2 to ~ —2.5 depending on the initial degree of 
mass segregation of the stars that have formed so far in the 
protocluster clump. 

Basu & Jones (2004) explored the time evolution of a 
log-normal distribution of the mass function of dense cores 
under the eflfect of accretion in a single zone model. They 
showed that by adopting a mass accretion rate formula of 
the type dM/dt oc M (or dM/dt oc M'^^'^) coupled to a 



nonuniform accretion timescale of the cores as a function 
of their mass, the mass function develops a power low tail 
in the intermediate to high-mass end while retaining its 
log-normal shape in the low mass regime. Bate & Bonnell 
(2005) presented an accretion-ejection scenario for the IMF 
based on the idea that all objects form with the same initial 
mass which is set by the choice of the opacity limit, and 
then accrete at a constant rate until they are ejected from 
their initial environment. The accretion rate in the Bate 
& Bonnell model are drawn from a lognormal distribution 
of a given mean accretion rate and an associated disper- 
sion. The ejection of the objects is described by a proba- 
bility exp{—t /reject), where reject is a given half life of the 
ejection process. Bate & Bonnell (2005) showed that their 
model can, with appropriate parameters, reproduce some of 
the standard Galactic field IMFs (i.e., notably the Miller & 
Scalo IMF). Myers (2009) showed that the evolution of the 
dense core mass function under the effect of accretion may 
also depend on the geometry of the environment in which 
the spherical star forming clump is embedded. In particular, 
he showed that the development of a flatter slope at the high 
mass end is strongly dependent on the environment density 
and weakly dependent on the environment dimensions. 

While our model shares a few conceptual similarities 
with the above described models, it also contains a number 
of noticeable different approaches. First, unlike the other 
models, our prescription for the time dependent accretion 
rate shape and normalization is parameter free since it is 
calibrated from numerical simulations. Another diflference 
to the previous models is the multi-zone (various location 
in the clump) and multi-epoch (cores injected over several 
epochs) aspect assigned to our treatment of the populations 
of PSCs and of stars and the co-evolution of two separate 
populations (only Bonnell et al. 2001b considered a multi- 
zone model but still not a multi-epoch model). The result of 
this dichotomy into PSCs and stars is that PSCs that have 
collapsed to form stars are subtracted from the local PSCs 
distributions. Also, unlike previous models, feedback from 
the newly formed stars in taken into account in order to pro- 
vide an estimate of the epoch at which gas is evacuated from 
the protocluster environment. The dichotomy into PSCs and 
stellar populations also allows us to investigate the effect of 
various core properties on the resulting IMF such as their 
timescale of contraction (parameter v) and their mass-peak 
density relation (parameter /i) which have not been taken 
into account by any previous model. Obviously the model 
presented in this paper can be further improved/modified 
in future works when more accurate theoretical and obser- 
vational constraints become available. Future improvement 
might be to adopt an accretion rate prescription that in- 
cludes the eflfects of the magnetic field which were not taken 
into account in the simulations of SK04. The eflfect of the 
magnetic field is, in addition to reduce the probability of 
core formation with increasing field strength (e.g., Dib et al. 
2007b), is also to strongly reduce the accretion rates onto 
the cores (Price & Bate 2008). Another modification would 
be to take into account the sub-fragmentation of the cores 
and build a single-star stellar IMF. As stated in § [6l we 
made in this work the assumption that a single core leads to 
a single star/star system and thus the derived IMFs can be 
compared to system IMFs rather than to single star IMFs. 
Choosing adapted accretion rates as a function of the degree 
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of magnetization of the clump may also allow us to simul- 
taneously select an adapted value of v for a given magneti- 
zation level. Further improvements may also include using 
more constrained values of the parameter k which describes 
the fraction of the wind energy that opposes gravity. Sim- 
ulations of multiple wind interactions are needed to further 
constrain its value. 



9 SUMMARY 

In this work, we use semi-analytical modeling to study the 
effects of gas accretion onto populations of dense gravita- 
tionally bound cores at different locations in a protocluster 
clump and the transition of the cores to form stars. Once 
the local contraction timescales for cores of a given mass 
are shorter than their accretion timescales, cores are turned 
into stars. We are thus able to follow the co-evolution of the 
dense core mass function (PSCMF) and of the IMF. Feed- 
back from massive stars in the form of stellar winds is taken 
into account. The energy from stellar winds that could dis- 
perse the gas from the protocluster clump is compared to its 
gravitational energy and once it becomes dominant, the pro- 
cess of gas accretion onto cores is terminated, and the star 
formation process is quenched, thus setting up the IMF. We 
rely on observations and direct numerical simulations in or- 
der to calibrate the properties of the protocluster clumps, 
their populations of cores that are prone to become stars, 
and the physical processes of gas accretion onto the cores 
and feedback from massive stars in the form of winds. Af- 
ter prescribing a clump mass, the main free parameters of 
this Accretion-Collapse-Feedback (ACF) model are a) the 
exponent of the clump density profile, 6, b) the contraction 
timescale of the cores in units of their free-fall timescale i/, 
and c) the exponent of the mass-concentration relation of 
the cores, fi. A fiducial case of the ACF model is discussed 
in detail and compared to the Orion star forming region. 
The fiducial case of the ACF model reproduces: 

a) The normalization and complex shape of the IMF of the 
Orion Nebula Cluster (ONC). Namely, the shallow slope of 
the IMF of the ONC in the mass range ~ [0.3 - 2.5] M©, 
a steeper slope in the mass range ~ [2.5 — 12] Mq, and a 
nearly flat tail at the high mass end. 

b) The age spread of 80 percent of the young stars in Orion. 
We find an age spread of ~ 2.2 x 10^, whereas Hillenbrand 
(1997) suggested that 80 percent of the stars in Orion are 
younger than 3 x 10^ years. 

c) The location of the most massive stars which are born in 
the region located between 2 — 4 Rco from the clump center, 
where Rco = 0.2 pc is the core radius of the clump in the 
fiducial model. 

d) Simultaneously with the IMF, the normalization and 
shape and of the mass function of the dense sub-millimeter 
cores in the star forming regions of Orion which displays a 
tail at the high mass end. 

In a second step we evaluate the effects of changing 
some of the models parameters on the resulting IMF. The 
results are the following: 

e) More concentrated cores with increasing mass (i.e., larger 
values of p in the models), result in steeper IMFs (Fig. llOf) . 
whereas smaller value of fi favor the formation of massive 



stars and the development of a tail-like structure at the high 
mass end. 

f) Longer contraction timescale of the cores (i.e., larger val- 
ues of V in the models) favor a longer and more efficient 
accretion onto the cores and thus the formation of massive 
cores and in the same time the depletion of lower mass cores. 
This results in IMFs which contain an increasing fraction of 
massive stars and no or very small fractions of smaller mass 
stars. On the other hand a decreasing contraction timescale 
of the cores favors the formation of a steeper IMF at the 
high mass end (Fig. I12p . 

g) Star formation in a protocluster clump with a radial pro- 
file which is close to an show little variations in the 
resulting IMFs. However, stronger deviations of the density 
profile of the clump result in significant variations both at 
the low and high mass end. A steeper radial profile of the 
clump results in an IMF which contains massive stars and a 
relatively depleted low mass end, whereas a shallower clump 
profile results in an IMF which is steeper at the high mass 
end, shallower at the low mass end, with a broad interme- 
diate mass plateau (Fig. I14p . 

h) Variations in the level of turbulence in the protoclus- 
ter clump (i.e., variations in the exponent of the velocity 
dispersion-size relation) essentially result in variations at the 
low mass end in the resulting IMFs. Higher levels of turbu- 
lence result in an IMF that contains larger fractions of low 
mass stars (Fig. I17|l . The shape of the IMF at the high mass 
end is unaffected by changing the turbulence levels in the 
clump. This indicate that the resulting IMF, in this mass 
regime, is not very sensitive to the exact shape of the mass 
function of dense cores that are being continuously gener- 
ated in the clump until star formation is quenched by stellar 
feedback. 

Overall, it is important to note that the IMF that re- 
sults from various combinations of the most relevant param- 
eters of the ACF model is usually not well described by a 
broken power law function or by a log-normal distribution. 
Instead, the resulting IMFs show complex shapes and gen- 
uine variations in all mass regimes that can extend beyond 
the statistical uncertainties. Provided there are variations in 
the structural and dynamical properties of real protoclusters 
clumps and in the populations of cores prone to star forma- 
tion within them, our model predicts genuine variations in 
the stellar IMF. 
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iVlcdlillit Wi Llitj Vdji lCLLJlt:l5 


Mc (P£ 


mass of the clump 




radius of the clump 




core radius of the clump 




velocity disperion of the gas in the clump 


PcO 


central density of the clump 


6 (!?;[£ 


exponent of the density profile of the clump at large radii 




gravitational energy of the protocluster clump 




exponent of the velocity dispersion-size relation in the clump 




exponent of the kinetic energy power spectrum in the clump 




free-fall timescale of the clump 




fraction of the mass of the clump that turn into dense cores per t//,c 


Core variables 


Meaning of the variables 


M 


mass of the core 


PpO 


central density of the core 


(IPlI 


exponent of the ppo — ^ relation of the core 


Rpo 


core radius of the core 


Rp 


radius of the core 


Pp 


average density of the core 


*ff 


tree-tall timescale ot the core 


icont ,p 


contraction timescale of the core 




ratio of the contraction timescale of the core to its free-fall time 


Accretion variables 


Meaning of the variables 




effective accretion rate of the core 


Stellar and Feedback variable 


Meaning of the variables 


« (pii 


fraction of the mass of the core that ends up in the stars 




masses of stars 




mass loss rates from stars 




time integrated energy from massive stars for M ^ 10 M© 


K (iPii 


fraction of E^ind converted into motions that oppose the gravity of the clump 


^k,wind 


time integrated energy from stellar winds that oppose the gravity of the clump 



The free parameters of the models are marked with (P), and the important ones are marked with (IP). Important 
parameters are the ones that can cause variation in the shape of the PSCMF and of the IMF whereas the other 
parameters will only affect the normalization (i.e., vertical shift) and the horizontal shift along the mass axis. Note 
that a is not a free parameter of the model as it is constrained by the observed — Rc relation by Saito et al. 
(2007). However, its value has been varied around the observed central values within the limits of the statistical 
uncertainties. 

^ In all the models, we chose Mc = 10* Mq which is a characteristic mass for protocluster clumps. To first order, 
changing will change the vertical normalization of the PSCMF and of the IMF but not their shapes. 

In all models, we chose R^o = 0.2 which is characteristic of protostellar clumps and of their stellar cluster 
progenitors (e.g., Hillenbrand & Hartmann 1998). 

In the observations, b is seen to vary between ~ 1.5 to 2.5. In our models b is varied between 1.4 and 2.2. 
^ The core formation efficiency or fraction of the mass of the clump that turn into dense cores per t//,ci Ec may vary 
from clump to clump depending on physical conditions such as the magnetic field strength in the clump ranging 
from 1 % or less, to a few tens of percent (e.g., Nakamura & Li 2008; Vazquez-Semadeni et al. 2005b ;Padoan & 
Nordlund 2009). In all the models Ec was set equal to 0.07. Note that similarly to Mc, the choice of Ec, as long 
as feedback is not considered, does not change the shapes of the PSCMF and of the IMF but only their vertical 
normalization. 

f In our models, the exponent of the ppo — ^ relation of the cores has been allowed to vary in the range 0-0.6 in 
agreement with the range determined from observations (e.g., Caselli & Myers 1995; Johnstone &: Bally 2006). 
9 The value of v is expected to vary between 1 and 10. Numerical and observational determinations of u suggest 
that it is in the range of a few (e.g.. Ward- Thompson et al. 2007; Galvan-Madrid et al. 2007). 

^ In all models, ^ has been fixed to the value of 0.1. Varying ^ does not change the shape of the resulting IMF but 
only shifts it horizontally along the mass axis. 

* In all models, k has been fixed to the value of 0.1. As explained in the text, all that is presently known is that 
K ^l. Simulations of multiple wind interactions arc needed to further constrain its value. 

Table 1. Main variables in the Accretion-Collapse-Feedback model. From top to bottom, the first panel describes the protocluster clump 
variables, the second panel the dense cores variables, the third panel the acrretion variable, and the lower panel the stellar and stellar 
feedback variables. 
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